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We compute (P(Qs) QCD corrections to the lepton invariant mass spectrum in the decay b —> ulvi, 
relevant for the determination of the CKM matrix element \Vub\. Our method can also be used to 
evaluate moments of the lepton energy distribution with an 0(0^} accuracy. The abelian part of 
our result gives the neutrino invariant mass spectrum in the muon decay and, upon integration, the 
O (q^) correction to the muon lifetime. 
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Determination of the CKM matrix elements from pre¬ 
cision studies in i?-physics is one of the main goals of 
experiments BaBar and Belle, under way at SLAG and 
KEK. These studies are expected to provide important 
insights into flavor physics, in particular to shed light 
on the origin of the CP-violation and, possibly, discover 
“New Physics.” 

The two CKM parameters directly accessible at B- 
factories, \Vcb\ and |lAb|, strongly differ in magnitude, 
with the former being about ten times larger than the 
latter. An accurate determination of \Vcb\ is much easier 
since the relevant decay rates are relatively large and the 
backgrounds are small. For |14b|, the theoretically favor¬ 
able methods are not feasible experimentally, whereas in¬ 
terpretation of clean experimental signatures suffers from 
large theoretical uncertainties. 

Extraction of \Vub\ from inclusive semileptonic decays 
of B mesons requires a suppression of the much larger 
contribution of 6 —> c transitions. In order to do so one 
has to impose cuts on various observables and several op¬ 
tions have been discussed in the literature. For example, 
one can select events with large energy of the charged 
lepton, which can be produced only in 6 ^ u decays, or 
require that the hadron invariant mass be smaller than 
the lightest charmed meson D. 

Unfortunately, such cuts are so severe that the rate 
of the remaining events cannot be predicted using the 
Heavy Quark Expansion. For example, imposing the cut 
on the electron energy induces a sensitivity of the decay 
rate to the H-meson light-cone wave function which is 
not very well known. It can in principle be extracted 
from measurements of the photon energy spectrum in 
b —> s"f. However, the relevant theoretical analysis has 


only been performed in the limit of an infinite b quark 
mass and the potentially sizable Aqcb/ mb corrections 
are not under control. It is desirable, therefore, to have 
an alternative combination of cuts which can remove the 
charm background, keep a significant fraction of 5 —> 
u events, and preserve the applicability of the standard 
Heavy Quark Expansion. 

Recently, a method fulfilling these requirements has 
been proposed by Bauer, Ligeti and Luke 01- Their 
idea consists in extracting |I4b| from inclusive semilep¬ 
tonic decays b ului by applying a cut on the invari¬ 
ant mass of the leptons q^. To eliminate the charm 
background, one requires q^ > q^ = {mB — mu)^ ~ 
11.6 GeV^. It turns out that this cut is mild enough to 
keep significant fraction of 6 —> u transitions and also 
the energy release is sufficiently large so that the stan¬ 
dard methods of the Heavy Quark Expansion in 1/mb 
can be applied with confidence. Of course, there are 
several sources of the theoretical uncertainties associated 
with this method, including in roughly equal measure the 
value of the b quark mass, the non-perturbative power 
corrections (of third order in the ratio of Aqcd and the 
characteristic momentum flow), and the two-loop pertur¬ 
bative QCD corrections 0. The calculation of this last 
effect is the main purpose of this Letter. 

The difficulty connected with such corrections is that 
they involve the q^ distribution, rather than the total de¬ 
cay rate. While two-loop corrections to charged particle 
decays are in general challenging (the first calculations for 
specific kinematic configurations or the total decay rates 
have been completed only recently, see e.g. i I i 0), 
two-loop corrections to the decay distributions have never 
been evaluated so far. 
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FIG. 1: Examples of diagrams whose cuts contribute to the 
semileptonic decay b uli/f. (a) abelian; (b) light or heavy 
quarks; (c) non-abelian. 


In the present calculation we take advantage of the 
fact that, for the experimentally interesting case, the in¬ 
variant mass of the leptons is large. We introduce an 
expansion parameter S = {m^ — q^)/ml. In b —>■ u stud¬ 
ies using cuts proposed in [Q ^ the maximal value of S 
is about 0.5 for q'^ — q^. Obviously, increasing q^ re¬ 
sults in a rapid decrease of (5, so that 5 can be considered 
as a small parameter in the region of interest q^ > q^. 
Therefore, by constructing an algorithm for expanding 
the relevant Feynman diagrams around 5 = 0 and com¬ 
puting several terms of such an expansion, we can derive 
the 0{a^) correction to the dilepton invariant mass spec¬ 
trum valid in the region of experimental interest. 

Examples of diagrams we have to consider in studying 
the semileptonic 5 —> n decay at O [al) are shown in 
Fig. 0. The optical theorem connects the imaginary part 
of such diagrams with contributions to the decay. We 
first integrate over the lepton and neutrino phase space, 
thereby reducing the problem to the decay b W*(q^)u, 
where W* is a virtual W boson with an invariant mass 
q^. In the limit S —>■ 0, q^ approaches mf. There¬ 
fore, due to phase space constraints, W* becomes static. 
The expansion in S is constructed by applying the Heavy 
Quark/Boson Expansion to the Feynman diagrams. The 
only unusual feature in our case is that the initial 6 -quark 
is on the mass shell. In the HQET limit, this leads to 
propagators of the type l/(2pk), whereas the W* boson 
is off-shell so that its propagator has the form l/(2pk+S). 

Since we are interested in the O (a^) corrections to the 
decay distributions, we have to consider the three-loop di¬ 
agrams of the self-energy type, like those shown in Fig. 0 
and extract their imaginary parts. Initially, there are 
two scales in the problem: using mb as a unit of energy, 
these scales can be expressed as O {1) and 0{5). We 
employ asymptotic expansions to indentify contributions 


arising from these widely separated scales. The region 
with all loop momenta of 0(1) does not contribute to 
the imaginary part since it is analytic (polynomial) in 
5. When some loop momenta are O (5) and others are 
0 (1), a three-loop diagram factorizes into a product of 
one- and/or two-loop diagrams and is easy to evaluate. 

The non-trivial part of the calculation is the HQET 
limit where all loop momenta are of 0(5). These di¬ 
agrams are similar to the three-loop HQET diagrams 
10 , 0, ^ but not identical with them, since some of 
the lines in the present case are on-shell. We have con¬ 
structed an algorithm based on recurrence relations and 
integration-by-parts identities [pd]| with which one can 
reduce any relevant three-loop diagram to a linear com¬ 
bination of a few master integrals. Four of these master 
integrals are new. We compute them in the Euclidean 
ijp' = —1), O = 4 — 2e dimensional space. 

Propagators occuring in the master integrals are de¬ 
noted by Di = kl, D 2 = A:|, D 3 = A:|, = {ki — 

D 5 = (fe - ksf, Dq = 2 pki, Dr = 2 pk 2 , Ds = 
2 pk 3 , Dq = 2 pki + 2 pk 2 , Diq = 2 pki + 2 pk 3 , Du = 
2pki + 2 pk 2 + 2 pk 3 . The four new results are 
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In the above formulas ((3 is the Riemann zeta function, 

00 

Cs = i] 1/*^, and Ce = + e). 

i—1 

Using recurrence relations to reduce all loop integrals 
to a combination of master integrals (these algebraic ma¬ 
nipulations are done with FORM |l^), we obtain the 
0{a'() correction to the dilepton invariant mass spectrum 
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(we use the pole mass nif, and the MS scheme for ag) 
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where L = \t\5 and Cf = 4/3, Ca = 3, and Tr = 1/2 
are the usual SU(3) color factors and Nr and Nr denote 
the number of light (m, = 0 ) and heavy {iriq = mb) 
quark species. We use the approximation rric = mb since 
for > Qq there is no phase space available for charm 
quark production. If needed, corrections for rric 7 ^ mt in 
virtual effects can easily be computed. 

For the coefficients Xa, -XIna, Xr, Xr we find 
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For brevity we have presented the results accurate up 
to the terms 0{S^). For the numerical analysis below we 
use terms up to 0((5®). 

We tested these results in several ways. We used a 
general covariant gauge and checked the cancellation of 
the gauge parameter. The result for Xr agrees with the 
numerical calculation in [p^ . A simple interpolating for¬ 
mula which we actually used for the comparison can be 
found in the appendix of Ref. . The agreement is very 
good, practically for all values of S. 

Further, we can extrapolate the results of the expan¬ 
sion by taking the limit <5 ^ 1 in which case our formu- 
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las should describe the decay of a massive quark into a 
massless quark and a massless W boson. In this limit, 
second order QCD corrections were computed for the top 
quark decay |^. We find that for the color structures 
-^NA,L,H the difference between the two results is better 
than 10%. The agreement is much worse for the abelian 
part Xa, where the difference can be as large as 50%. 
This demonstrates that the seven terms of the expansion 
are insufficient for the abelian part to converge in the 
limit (5 —> 1. 

However, because of the SU(3) color factors, the con¬ 
tribution of the abelian part is suppressed and we can 
reliably derive the 0(ag) correction to top quark decay 
from our formulas. Taking Xl = 5 and Xff = 1, we 
find X 2/2 « —16.4, whereas the central values of the co¬ 
efficients in give —16.7. An even better agreement 
is obtained for S = 6w = 1 — M^/m^ ci 0.79. At this 
point, corresponding to physical values of the W boson 
and top quark masses, the width of < —> bW was evalu¬ 
ated in iQ . We have perfect agreement with the central 
value of the second order correction, X2{Sw)/‘2 = —15.6, 
given in eq. (28) of that paper. 



FIG. 2: O (os) correction to the decay width b ulvi, 
X 2 {S) (defined in eq. (|^), as a function of 5 = 1 — Iml (for 
Nl = 3, Nh ~ 2). 

As the final check one can integrate Eq. i) over (5, 
obtaining the total decay rate b —>■ ului, for which the 
second order QCD corrections are known |^. Taking 
A'l = 4 and Nh = 1 and integrating over 6 we ob¬ 
tain /p d(5A2((5) = -21.24, in excellent agreement with 
—21.296, given in Ref. 0. 

Integrating the abelian contribution Xa we can com¬ 
pute the two-photon corrections to the muon lifetime. 
We find dSX aIS) ~ 3.1, where the 13% discrepancy 
with the exact value in eq. (9) of is due to poor con¬ 
vergence of our series for large 6. However, if we assume 
that the convergence is good up to (5 ~ 0.65 and extrapo¬ 
late for larger S using Aa( 1) = 7.0(4) ||^, we reproduce 
the muon lifetime correction within 3%. 

For S ^ 1/2, relevant for the extraction of \Vub\, the 
series converge very well and accurately approximate all 
color components of the O (a^) correction. 


The full O (ttg) correction to the quark decay width, 
X 2 {S), is plotted in Fig. |[ Even at the end point <5=1, 
our estimate for X 2 agrees with our result for the top 
decay to better than 3%. 

To show the impact of the computed corrections on 
dilepton invariant mass distribution, we separate the 
BLM and non-BLM corrections since the former have 
already been studied in the literature. We define the 
BLM and non-BLM corrections as 

^BLM ^ -SCfPqXl, = A 2 - Xf^^, (5) 

where/3o = 11(7^/12—TRAfi/d denotes the beta-function 
coefficient in a theory with three massless quark flavors, 
appropriate for the range of used for the I4b extrac¬ 
tion. 



FIG. 3: The non-BLM corrections, Ar“®^'^/(6<5^ - 4<5^) (for 
Nl = 3, Nh = 2). 

The value of the BLM corrections is known to be 
strongly correlated with the scale of the coupling con¬ 
stant used in the one-loop result and also with the quark 
mass used in the formula for the decay rate A discus¬ 
sion of these issues can be found in the literature @ and 
we will not consider them here. On the contrary, the 
non-BLM corrections are new. Their dependence on S 
is shown in Fig. ^ where the ratio of the non-BLM cor¬ 
rections and the tree level decay rate 65^ — 4(5^ is plot¬ 
ted. For realistic values of the strong coupling constant, 
as = 0.2 —0.3, the non-BLM corrections are about 5% in 
the range of S relevant for the \Vub\ extraction from the 
dilepton invariant mass spectrum. 

The technique described in this Letter might open a 
way to reliable estimates of the O (a^) corrections to 
more complicated observables. For example, a simple 
modification allows one to calculate the moments of the 
charged lepton energy spectrum for a fixed value of the 
dilepton invariant mass. 

Recently, combined cuts on both dilepton and hadron 
invariant masses were advocated for the \Vub\ determi¬ 
nation iQ . It has been argued that in this approach one 
can keep the theoretical uncertainties under control while 
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retaining a larger data sample of the b u transitions. 
Since the calculation reported here has been performed 
without any restriction on the hadronic invariant mass, 
our results for the QCD corrections are not applicable 
in this case. However, a sufficiently large number of 
moments should contain enough information about the 


spectrum to determine the effect of the cut. 
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